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The  critical  level  for  ozone,  above 
which  it  has  a detectable  effect  on  bio- 
logical targets,  is  potentially  to  be  set 
by  the  United  Nations  Economic  Com- 
mission for  Europe  at  300  nL-h/L 
hours  per  annum  over  40  nL/L.  It  is 
therefore  important  to  determine  the 
aggregate  exceedance  over  40  nL/L 
throughout  the  United  Kingdom.  Over 
most  of  the  UK,  ozone  concentrations 
are  unknown  so  we  rely  on  our  under- 
standing of  the  atmospheric  processes 
and  on  the  statistical  properties  of 
ozone  concentrations  to  interpolate 
between  monitoring  sites.  This  paper 
describes  the  application  of  statistical 
models  derived  for  storm  severity  data 
to  the  ozone  data  for  the  United 


Kingdom.  Aggregate  excess  distributions 
were  fitted  to  data  from  all  rural  moni- 
toring sites  using  a Weibull  model  with 
a 40  nL/L  threshold.  At  this  threshold 
the  scale  parameter  has  a spatial  inter- 
pretation, but,  with  higher  thresholds, 
there  were  problems  with  missing  data 
and  small  scale  spatial  effects  were  not 
detected.  The  approach  appears  suc- 
cessful for  all  except  very  large  aggre- 
gate exceedances  which  deviate  from 
the  Weibull  predictions. 
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1.  Introduction 

The  major  public  concern  with  ozone,  O3,  in 
Europe  has  focused  recently  on  the  existence  of 
“ozone  holes”  in  the  stratosphere  caused  by  the  de- 
pletion of  ozone  as  a consequence  of  chlorofluoro- 
carbon emissions.  Ozone  is  also  present  in  the 
troposphere  and  in  the  planetary  boundary  layer  at 
concentrations,  i.e.,  volume  fractions  between  10 
and  200  nL/L  (i.e.,  parts  per  billion,  ppb  = 10”'').  In 
the  second  half  of  the  last  century  European  mean 
concentrations  ranged  between  10  ppb  and  15  ppb 


[8].  Current  mean  concentrations  are  about  twice 
these  values  and  ozone  episodes  with  peak  concen- 
trations between  100  ppb  and  200  ppb  occur,  a level 
known  to  cause  damage  to  many  plant  species. 
Episodes  happen  if  the  precursor  gases  for  photo- 
chemical ozone  production  (oxides  of  nitrogen,  NO 
and  NO2,  and  volatile  organic  compounds,  VOCs) 
are  present  in  suitable  meteorological  conditions 
for  the  chemical  reactions  to  occur  (ideally  hot  sum- 
mer days  with  clear  skies  and  low  wind  speeds). 
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The  description  of  spatial  patterns  in  exposure  of 
vegetation  to  ozone  over  Europe  has  been  ham- 
pered by  the  limited  availability  of  monitoring  data, 
the  very  large  spatial  variability  in  ozone  concentra- 
tions and  a poor  understanding  of  the  underlying 
mechanisms  regulating  the  ozone  exposure  of  ter- 
restrial ecosystems.  Defining  a threshold  for  phyto- 
toxicity is  not  simple.  However,  at  60  ppb  of  ozone 
there  is  little  doubt  that  there  is  a clear  contribution 
from  photochemical  production  in  polluted  air  and 
a map  of  hours  over  60  ppb  for  Europe  (Fig.  1)  is  a 
guide  to  some  broad  trends  [4].  In  a large  area 
north  of  the  Alps,  covering  most  of  Germany  and 
parts  of  neighbouring  countries,  200  hours  per  year 
above  60  ppb  is  common.  North  and  west  of  this 
area  the  annual  duration  of  exposure  declines  but 
to  the  east  there  is  so  little  information  available 
that  mapping  is  uncertain.  The  Mediterranean  zone 
of  high  ozone  exposure  reflects  recent  work 
showing  that  ozone  episodes  are  common  events  in 
this  region  but  the  levels  have  not  yet  been  well 
quantified. 


Fig.  1.  Hours  when  ozone  exceeds  60  ppb. 


Although  the  meteorological  conditions  leading 
to  ozone  episodes  are  similar  in  a general  sense  at 
all  sites,  the  climates  of  northern  and  southern  Eu- 
rope lead  to  very  different  patterns  of  events.  In 
northern  Europe,  typical  episodes  occur  when  a sta- 
tionary spring  or  summer  anticyclone  provides  the 
conditions  for  ozone  production  from  the  emitted 
precursor  gases  to  add  appreciably  to  the  back- 
ground concentration  of  about  30  ppb.  Typical  pro- 
duction rates  give  net  increases  of  10  ppb  to  20  ppb 


per  day  and  a succession  of  8 to  10  such  days  leads 
to  peak  concentrations  of  150  ppb  to  200  ppb.  Often 
in  northern  Scandinavia,  Britain,  Ireland  and  west- 
ern France  the  ideal  meteorology  exists  but  in  the 
absence  of  upwind  precursors.  In  Germany  and 
central  Europe  almost  all  wind  directions  provide 
the  precursors  and  hence  the  NW-SE  gradient  in 
episodes.  In  southern  Europe,  the  meteorological 
conditions  are  more  stable  and  episodes  can  occur 
every  day  for  long  periods.  However  the  effects  of 
both  sea  breezes  and  the  development  of  intense 
thermal  low  pressure  areas  on  the  air  circulation 
causes  very  variable  patterns  of  ozone  exposure. 

Superimposed  on  this  two  dimensional  surface 
there  is  a daily  cycle  in  ozone  concentration  which 
is  a very  important  and  variable  feature.  At  low  al- 
titude inland  sites  a marked  diurnal  variation  (of 
the  order  of  30  ppb)  is  observed  but  at  high  eleva- 
tion the  amplitude  of  the  diurnal  cycle  gradually  re- 
duces to  less  than  5 ppb  at  mountain  tops.  Figure  2 
shows  data  for  1 day  at  both  Great  Dun  Fell  (847  m 
above  sea  level)  and  Wharleycroft  (206  m above  sea 
level),  two  sites  which  are  less  than  10  km  apart. 
Hill  tops  are  generally  windy  sites  at  which  the  ter- 
restrial surfaces  are  well  connected  to  the  free  tro- 
posphere and  where  the  downward  supply  of  ozone 
to  the  surface  exceeds  the  rate  of  deposition.  At  low 
level  sites  the  thermal  stratification  of  the  atmo- 
sphere with  the  development  of  a nocturnal  inver- 
sion restricts  the  supply  of  ozone  from  above  during 
the  night  and  morning.  In  these  conditions  both  de- 
position to  the  surface  and  the  nocturnal  atmo- 
spheric chemical  titration  of  ozone  with  nitric  oxide 
causes  the  surface  concentrations  of  ozone  to  de- 
cline, potentially  to  negligible  levels.  At  coastal  sites 
the  effects  of  land  and  sea  breezes  strongly  modify 
the  ozone  exposure  of  the  ground. 


Fig.  2.  Typical  altitude  effect  on  the  diurnal  cycles  of  ozone 
concentration. 


354 


Volume  99,  Number  4,  July-August  1994 

Journal  of  Research  of  the  National  Institute  of  Standards  and  Technology 


The  main  concern  with  rural  ground  level  ozone 
concentrations  is  the  damage  which  can  be  caused 
to  plants  and  to  human  and  animal  health. 

For  vegetation,  some  sensitive  species  show  visi- 
ble or  physiological  effects  following  exposure  to  40 
ppb  or  50  ppb  [3].  However,  the  effect  of  exposure 
can  be  modified  by  the  presence  of  other 
atmospheric  pollutants  and,  since  ozone  causes 
damage  to  vegetation  through  stomatal  uptake,  by 
nutritional  status,  light,  temperature  and  humidity. 
There  is  genetic  variability  in  the  ozone  response 
within  species  as  well  as  between  species  and,  al- 
though considerable  attention  has  focused  on  crops 
and  forests,  little  is  known  about  the  impact  on 
semi-natural  vegetation.  Timing  of  the  exposure 
within  the  life  of  the  plant  can  be  important  as  can 
be  the  time  for  recovery  between  exposures  [6].  The 
United  Nations  Economic  Commission  for  Europe 
is  considering  a tentative  proposed  critical  level  of 
300  ppb  • h above  40  ppb  during  daylight  hours  for 
the  growing  season  of  the  vegetation.  A critical  level 
is  defined  as  one  below  which  ozone  has  no  de- 
tectable effect.  However,  there  are  a number  of  out- 
standing issues  which  it  is  hoped  to  resolve  by  the 
end  of  1993  and  the  adopted  critical  level  may  well 
be  different.  The  proposed  critical  level  would 
probably  be  exceeded  in  most  of  Europe  at  present. 

Concern  for  human  health  in  the  UK  at  the  cur- 
rent levels  of  ozone  is  growing  but  better  assess- 
ments of  population  and  individual  exposure  are 
thought  necessary  [6].  This  aspect  may  in  time  be 
the  main  argument  for  emission  controls  of  the  ma- 
jor precursor  gases. 

It  is  important  to  differentiate  between  the  dose 
which  a plant  or  human  receives,  that  is  incorpo- 
rated into  the  individual’s  system  by  some  method, 
and  its  exposure,  that  is  the  level  in  the  atmosphere 
around  the  individual.  In  this  paper  current  meth- 
ods for  determining  plant  exposure  in  the  UK  are 
described  and  then  the  potential  application  of  ex- 
treme value  theory  is  explored. 

2.  Ozone  Exposure  Maps  of  the  UK 

Ozone  exposure  has  recently  been  mapped  for 
the  UK  at  three  concentration  thresholds:  40  ppb, 
60  ppb  and  90  ppb  [5].  There  were  about  17  rural  or 
semi-rural  monitoring  stations  between  1987  and 
1991  which  recorded  hourly  mean  concentrations  in 
Britain  and  Ireland  (Fig.  3).  As  the  differences  be- 
tween sites  which  were  geographically  close  was  as 
large  as  the  differences  between  geographically  dis- 
tant sites,  a straight  spatial  interpolation  between 
sites  gave  a map  similar  to  that  in  Fig.  4. 


In  the  summer  months,  taken  as  April  to  Septem- 
ber, during  the  part  of  the  day  when  the  atmo- 
spheric boundary  layer  was  well  mixed  by 
turbulence,  ozone  concentrations  at  neighbouring 
sites  were  very  similar.  For  each  threshold,  an  em- 
pirical relationship  was  derived  between  the  hours 
over  the  threshold  for  the  whole  day  and  the  hours 
over  the  threshold  for  the  well  mixed  period,  taken 
to  be  1200  to  1800  GMT.  For  the  threshold  at  60 
ppb,  the  relationship  was 

/76,.  = ( 1.3 + 0.0021  2)  tfto,  (1) 

where  hm  was  the  total  hours  over  60  ppb,  twi  was  the 
hours  over  60  ppb  between  1200  and  1800  GMT 
and  2 was  the  altitude  of  the  location  in  meters.  This 
relationship  was  applied  to  the  spatial  interpolation 
of  hours  over  60  ppb  for  1200  to  1800  GMT  to 
provide  a map  (Fig.  5)  with  clear  topographical  in- 
fluence. The  coastal  effect,  which  can  extend  for  5 
km  to  20  km  inland  depending  on  meteorological 
conditions,  was  ignored;  typically  coastal  ratios 
were  around  2 rather  than  1.3.  The  maps  are  only 
for  the  summer  months,  April  to  September,  but 
ozone  levels  very  rarely  exceed  40  ppb  during  the 
remainder  of  the  year. 

This  approach  emphasizes  the  spatial  variability 
of  ozone  exposure  within  small  areas.  The  relation- 
ships for  the  different  thresholds  are  empirical  and 
must  be  recalculated  for  each  threshold  and  they  do 
not  provide  a general  description  of  high  concentra- 
tion events.  Direct  estimates  of  exposure  in  terms  of 
a dose  measurement  like  ppb. hours  are  not  avail- 
able although  a minimum  estimate  could  be  made. 
If  the  decision  were  made  to  set  different  windows 
to  match  the  growing  seasons  of  different  vegeta- 
tion types,  the  whole  procedure  could  be  difficult  to 
implement. 

3.  Modelling  Aggregate  Excess 

In  work  on  flood  levels  for  the  River  Thames,  An- 
derson and  Dancy  [2]  modelled  the  aggregate  ex- 
cess, that  is  the  sum  of  the  exceedances  over  a 
threshold,  within  a cluster  using  a Weibull  distribu- 
tion. There  are  similarities  between  ozone  data  and 
flood  level  data.  A Pareto  distribution  has  been 
shown  to  predict  the  peak  excesses  of  ozone  con- 
centrations at  a rural  site  using  a threshold  of  40 
ppb  [7].  There  is  a seasonal  component  in  the  data, 
since  high  concentrations  rarely  occur  over  the  win- 
ter period,  but  this  has  not  been  modelled  at 
present.  Also  ignored  was  the  probable  increase  in 
mean  values  of  ozone  concentration  over  the  time 
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Fig.  3.  Locations  of  15  monitoring  sites  on  the  UK  mainland 
(2  sites,  Lough  Navar  and  Mace  Head,  are  on  Ireland). 


Fig.  4.  Interpolated  map  of  the  number  of  hours  when  ozone 
exceeds  60  ppb. 


period  of  data  collection,  as  this  increase  was  small 
compared  to  both  the  diurnal  fluctuations  and  the 
accuracy  of  the  recording  methods. 

Anderson  [1]  has  looked  at  ozone  data  for  one 
site,  Stevenage,  for  a longer  time  period.  There  was 
evidence  of  nonstationarity  in  that  data  set  and  he 
shows  that  there  is  a need  for  temperature,  or  some 
similar  measure,  as  a covariate.  This  problem  is  still 
under  investigation  but  for  the  time  period  consid- 
ered in  this  paper,  1986  to  1991,  no  covariate  has 
been  used.  Anderson  also  derived  a method  of 
extrapolating  to  higher  thresholds  than  those  used 


in  fitting  the  models,  a very  useful  tool  for  deter- 
mining exposures  to  plants  with  different  sensitivi- 
ties to  ozone. 

The  data  for  the  17  sites  have  been  fitted  using  a 
single  threshold  of  40  nL/L  and  a Weibull  distribu- 
tion for  the  aggregate  excess.  The  data  were  declus- 
tered  using  a minimum  time  separation  of  48  hours. 
There  were  about  100  clusters  for  the  sites  with 
relatively  complete  data  sets.  Some  sites  were  not 
operational  in  the  earlier  period  of  collection  and 
one  site  had  only  3 years  data. 
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Fig.  5.  Altitude  adjusted  interpolated  map  of  the  number  of 
hours  when  ozone  exceeds  60  ppb. 

The  two  parameter  Weibull  model 

P{S  >s)  — Q\p{  — a'S*^)  (2) 

was  fitted.  The  shape  parameter,  9,  varied  between 
0.4  and  0.6  for  all  sites.  When  6 was  constrained  to 
the  value  0.5,  there  were  only  slight  increases  in  the 
values  of  the  likelihood  function.  The  spatial  varia- 
tion was  therefore  explored  using  only  the  scale 
parameter  a. 

The  Q-0  plots  showed  that,  as  expected,  the  fit  of 
the  Weibull  model  varied  from  site  to  site.  There 
were  some  very  straight  line  plots  but  there  were 


also  shapes  typically  illustrated  by  the  plot  for  Lady- 
bower  (Fig.  6).  Most  of  the  data  were  on  a reason- 
ably straight  line  but  the  Weibull  distribution 
underpredicted  a few  data  points,  usually  no  more 
than  five,  at  the  higher  values. 

The  shape  parameter,  a,  from  the  fitted  Weibull 
model  was  clearly  related  to  a SE-NW  trend  across 
the  country.  To  investigate  this  further,  the  a values 
were  regressed  on  other  available  data. 


Fig.  6.  Q-Q  plot  for  Weibull  model  (40  ppb  threshold,  48  hour 
separation)  fitted  to  data  from  Ladybower. 

The  sites  were  referenced  to  a line  from  Lulling- 
ton  Heath,  a site  on  the  south-east  coast  of 
England,  to  Strathvaigh  in  the  north  of  Scotland 
using  two  variables,  rimoi,  the  distance  north-west 
along  the  transect,  and  Cmod,  the  perpendicular  dis- 
tance from  the  transect  with  positive  values  being  to 
the  east.  The  actual  distances  were  divided  by  the 
length  of  the  transect  to  give  manageable  numerical 
values.  The  relationship  between  a and  «mod  was 
non-linear  and  the  simplest  function  of  rimod  which 
fitted  well  was  Umo^.  A linear  function  of  Cmod  im- 
proved the  fit.  One  site,  Bottesford,  had  a high 
residual.  This  site  has  peculiar  local  features  which 
can  give  it  the  characteristics  of  an  urban  site  and 
was  removed  from  the  data.  The  subsequent  regres- 
sion equation 

a = 0.091  + 0.083  n„,od  - 0.023  Cn,od  (3) 
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explained  about  85%  of  the  variation  and  produced 
an  acceptable  residual  pattern.  The  2 remote  sites, 
Strathvaigh  and  Mace  Head  (on  the  west  coast  of 
Ireland)  had  the  most  influence  on  the  fit. 

When  a threshold  of  60  ppb  was  used,  there  were 
problems  in  fitting  the  model  to  the  data.  Each  site 
had  only  about  four  independent  clusters  per  year 
and  fitting  a Weibull  model  was  very  difficult  since 
the  likelihood  surface  was  quite  flat.  Results  were 
obtained  by  assuming  6 was  the  same  for  a 60  ppb 
threshold  as  for  a 40  ppb  threshold.  However,  it 
became  apparent  that  there  were  potentially  two 
difference  sources  for  the  higher  exceedances  and 
that  separation  of  these  sources  was  critical  with  the 
increased  threshold  level. 

4.  Discussion 

The  results  of  these  fits  are  encouraging  although 
a number  of  problems  have  occurred.  The  Weibull 
distribution  with  a threshold  of  40  ppb  and  9 fixed 
at  0.5  gives  an  interpretable  underlying  pattern  for 
the  whole  country.  The  SE-NW  gradient  would  be 
expected.  Areas  to  the  east  of  the  chosen  transect 
are  more  influenced  by  air  masses  from  continental 
Europe  and  would  be  expected  to  have  more  ozone 
episodes.  The  lack  of  detection  of  an  altitude  effect 
at  this  threshold  is  not  entirely  surprising  as  high  al- 
titude sites  can  have  mean  ozone  concentrations 
quite  close  to  this  threshold.  However  an  altitude 
effect  would  be  expected  at  a higher  threshold. 

At  the  60  ppb  threshold  two  main  problems  oc- 
cur. The  first,  and  possibly  the  most  important,  is 
lack  of  data.  The  declustering  algorithm  which  has 
been  used  takes  the  rather  simple  approach  of  re- 
moving clusters  with  missing  data.  When  monitor- 
ing stations  are  running  continuously,  usually 
recording  several  times  per  minute,  there  are  a 
whole  series  of  glitches  which  can  occur  in  the  data 
for  reasons  wholly  unconnected  with  the  concentra- 
tion values.  In  particular  there  may  well  be  a series 
of  instrumentation  tests  which  usually  occur  during 
the  working  day  and  often  at  least  once  per  week. 
Data  capture  rates  of  over  90  % on  hourly  values 
are  regarded  as  good  but  not  all  sites  on  the  net- 
work are  achieving  these  rates.  Therefore,  careful 
decisions  on  the  treatment  of  missing  data  are  likely 
to  give  more  information  for  analysis. 

The  second  problem  is  one  of  determining 
whether  there  are  two  distinct  distributions  re- 
quired to  model  threshold  exceedances  or  whether 
the  Weibull  model  is  the  wrong  approach.  If  a mete- 


orological covariate  is  introduced  [I],  it  is  not  clear 
where  it  should  be  measured.  Rapid  ozone  produc- 
tion can  be  occurring  50  km  or  100  km  downwind  in 
good  sunny  conditions  but  the  monitoring  site  may 
be  sitting  in  quite  a different  climate.  Clearly,  fur- 
ther investigation  of  the  peak  and  close  to  peak  val- 
ues will  be  required. 

Although  large  exposure  to  ozone  can  be  accu- 
mulated by  a plant  at  concentrations  over  100  ppb, 
these  are  relatively  rare  occurrences  in  the  UK  and 
are  often,  if  not  always,  associated  with  very  dry 
conditions.  How  much  of  the  ozone  will  enter  the 
plant’s  system,  given  that  the  plant  is  probably 
under  considerable  water  stress  by  the  afternoon 
period,  is  not  clear.  Even  if  the  models  do  not  per- 
form very  well  at  the  highest  exceedances,  if  they 
can  perform  reasonably  well  for  the  remainder  of 
the  exceedances  they  could  be  of  considerable  ben- 
efit when  critical  levels  for  vegetation  are  consid- 
ered. For  human  health  problems,  of  course,  a 
different  perspective  is  required. 

This  approach,  when  combined  with  a model  of 
time  between  clusters,  has  the  potential  of  produc- 
ing valuable  information  for  the  assessment  and 
mapping  of  critical  levels  for  vegetation.  However, 
some  further  progress  is  required  with  models  for 
the  60  ppb  threshold  and  with  identification  of  local 
scale  variability  in  concentration  levels. 
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